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Abstract 

Clustering of cloud microphysical conditions, such as 
liquid water content (LWC) and drop size, can affect the 
rate and shape of ice accretion and the airworthiness of 
aircraft. Clustering may also degrade the accuracy of cloud 
LWC measurements from radars and microwave radiom- 
eters being developed by the government for remotely 
mapping icing conditions ahead of aircraft in flight. This 
paper evaluates spatial clustering of LWC in icing clouds 
using measurements collected during NASA research 
flights in the Great Lakes region. We used graphical and 
analytical approaches to describe clustering. The analyti- 
cal approach involves determining the “average size" of 
clusters and computing a clustering intensity parameter. 
We analyzed flight data composed of 1-s-frequency LWC 
measurements for 12 periods ranging from 17.4 minutes 
(73 km) to 45.3 minutes (190 km) in duration. Graphi- 
cally some flight segments showed evidence of consis- 
tency with regard to clustering patterns. Cluster intensity 
varied from 0.06, indicating little clustering, to a high of 
2.42. Cluster lengths ranged from 0.1 minutes (0.6 km) 
to 4.1 minutes (17.3 km). Additional analyses will allow 
us to determine if clustering climatologies can be devel- 
oped to characterize cluster conditions by region, time 
period, or weather condition. 

Introduction 

The National Aeronautics and Space Administration 
(NASA), the Federal Aviation Administration (FAA) Avia- 
tion Weather Research Program, and the Army Engineer 
Research and Development Center/Cold Regions Re- 
search and Engineering Laboratory (ERDC/CRREL) are 
assessing and developing radar and microwave radiom- 
eter technologies for remotely mapping icing conditions 
ahead of aircraft in flight. These systems are intended to 
increase aviation safety and utility by enabling aircraft to 
avoid potentially hazardous inflight icing conditions. Es- 
tablishing specifications for such systems requires con- 
sideration of limitations imposed by the operational en- 
vironment, by the radar and microwave radiometer tech- 


nologies, and by meteorological conditions. Understand- 
ing both the range and the variability of icing cloud 
microphysical parameters is critical for assessing radar 
and radiometer remote sensing performance and for es- 
tablishing specifications for the instruments. Examina- 
tion of the variability of icing cloud liquid water content 
(LWC) with clustering theory is the thrust of this paper. 

Background 

Clustering of cloud microphysical properties refers to 
clumping, or patchiness, in the spatial distribution of cloud 
microphysical parameters. Clustering implies that con- 
secutive values of LWC will be correlated over some dis- 
tance, defined as the spatial coherence length. If the mean 
of the entire data series is used to define clustering, then 
within a cluster, all the LWC values will tend to be greater 
than the mean or all values will tend to be less than the 
mean. Unlike a Poisson spatial distribution, where the 
parameters in the series are independent, in a clustered 
data series the values are not independent. 

Clustering of cloud microphysical properties may af- 
fect the design and accuracy of airborne or ground-based 
remote sensing systems designed to map icing conditions 
ahead of an aircraft. For example, experiments with mul- 
tiple-wavelength radar to retrieve cloud LWC 1-3 have 
shown that clustering can degrade the quality of LWC 
retrievals and can even cause reporting of negative LWCs. 
If explicit detection of LWC clusters is necessary to en- 
hance safety, shorter radar range gates may be required. 
Regardless, it will be necessary for radars to operate suc- 
cessfully in the presence of clusters. Within clustering 
conditions it may be desirable to also extend radar range 
to provide more advanced warning. Clustering may also 
affect the range, especially of shorter wavelength radars, 
because of the likely attenuation problems within clusters. 

It is not clear how radiometers will be affected by clus- 
tering of cloud microphysical parameters. Radiometers 
measure integrated brightness temperature and, being 
passive instruments, are not range-gated. As a result, ra- 
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diometers probably cannot detect clustering but also may 
not be affected by clustering. However, the ultimate range 
that a radiometer can “see” into clouds is a function of 
cloud attenuation caused by LWC, which could be af- 
fected by clustering. 

LWC and droplet size clustering can affect the rate 
and, as importantly, the shape of ice that forms on airfoils 
of fixed-wing aircraft. Aircraft traversing portions of cloud 
with small droplets and lower LWC in cold temperatures 
generally accrete rime ice. However, if extended clusters 
of higher LWC are intercepted, then clear ice may form 
with ice shapes, such as horns, that reduce lift, increase 
drag, and decrease control authority. This transition from 
rime to clear, or the reverse if LWC decreases, is a func- 
tion of the Schumann-Ludlam limit, which describes 
changes in LWC, drop size, and temperature magnitudes 
required to change the type, the density, and thus the shape 
of ice that forms on an airfoil. 4-6 Assessment of aircraft 
performance using models or icing wind tunnels simulat- 
ing non-clustered clouds may suggest that aircraft per- 
formance is acceptable when, in actuality, conditions may 
be hazardous if the same average cloud LWC were clus- 
tered. 

Though important, the clustering of cloud microphysi- 
cal parameters has not received as much attention as have 
other elements of cloud microphysics. Cooper et al. 7 sum- 
marized exceedance of LWC thresholds by season in Cali- 
fornia, Montana, Florida, Utah, Illinois, Kansas, and the 
Great Lakes region. Though not a true cluster analysis, 
the work indicates encounter distances and thus the maxi- 
mum distance over which certain LWC magnitudes were 
observed. Cober et al. 8 , using a similar procedure, de- 
scribe the patchiness of supercooled liquid water encoun- 
tered in the Canadian Atlantic Storms Project (CASP) 
using histograms of encounter number versus duration 
for patches of LWC greater than 0.025 g m -3 . 

A considerable amount of research has addressed the 
fluctuations of droplet size and LWC at small scales within 
clouds to assess turbulence, cloud evolution, and the ra- 
diative properties of clouds 9-17 . However, at the cloud 
scale of interest to us, the most appropriate analytical 
methods are presented by Jameson and Kostinski, in a 
series of papers beginning in 1997 4 13 . They describe clus- 
tering of drop sizes in rainfall and clouds using coher- 


ence length and a cluster parameter. The most recent pa- 
per 4 describes application of their techniques to LWC in 
synthetically generated icing clouds. This series of pa- 
pers is most applicable to our clustering problem, and it 
strongly influenced the work presented in this paper. 

Methodology 

We selected two techniques for assessing clustering — 
one graphical and one analytical. The intent of the graphi- 
cal technique is to allow visualization of the variability 
of the LWC along the flight path for specified LWC val- 
ues. The intent of the analytical method is to allow quan- 
tification of clustering and, perhaps, development of clus- 
ter pattern climatologies. These climatologies may be later 
inverted to develop correlated time series of LWC char- 
acteristics of specific cloud dynamics, such as overrun- 
ning along a warm front. 

Graphical Technique 

Graphically we used an incremented threshold ap- 
proach to show the frequency and percentage of time dur- 
ing a flight transect through a cloud when LWC exceeds 
specific threshold values. Higher frequencies of LWC 
crossings from below the threshold to above the thresh- 
old suggest more LWC variability at the specified thresh- 
old value along the flight track. If the threshold chosen 
coincides with the Schumann-Ludlam limit, then the fre- 
quency and duration of potentially dangerous icing con- 
ditions are identified. Jameson and Kostinski 4 suggest that 
more dangerous icing conditions can occur when an air- 
craft is flying in and out of areas that exceed the 
Schumann-Ludlam limit than when the aircraft is flying 
in an area that is Poissonian and does not exceed the 
Schumann-Ludlam limit for extended periods. 

For example. Figure 1 shows a small portion of 
a synthetically generated flight series with LWC values 
occurring at 1-s intervals generated from a random series 
with no clustering. Incremented LWC thresholds at 
0.02-g m 3 intervals from the same data yielded a plot 
(Fig. 2) with a linear decrease in percentage of time greater 
than the threshold as the threshold increases. The number 
of crossings increases from zero at a threshold smaller 
than the minimum LWC, to a maximum at the mean LWC 
for the data series, to zero again when the threshold be- 



1-s Measurements 


Figure 1. Synthetic random series with no clustering. The white line 
represents the mean of the data series. 
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Liquid Water Thresholds (g m ' 3 ) 


Figure 2. Cumulative distance and number of crossings of LWC 
greater than thresholds for nonclustered, Poissonian random 
series. 


comes larger than the maximum LWC in the series. The 
linear decrease of percentage of time exceeded by the 
thresholds as the threshold value increases, and the sym- 
metry of the number of crossings about the mean LWC, 
is a result of the randomness and lack of clustering of the 
series (Fig. 2). As illustrated later, flight segments of ac- 
tual cloud LWC clusters show different percentages of 
time greater than the threshold and different crossing fre- 
quency relationships. 


Analytical Technique 

We utilize the methodology of Jameson and Kostinski 4 
to analytically describe LWC clustering. If the cloud LWC 
is not clustered, it is Poissonian. The first step in the ana- 


lytical approach involves determining if clustering occurs 
in the data series. Kostinski and Jameson 13 indicate that 


a two-point autocorrelation function, as given in equation 
1, can be used to determine the presence of clustering: 



where r| (/) = two-point autocorrelation function 
(j, = mean of the series 

k( l )and k(0) = values at the reference location (time) and 
distance (time), respectively, at a distance 
(seconds) l from the reference. 


If clustering does not occur, that is, all values in the series 
are statistically independent, or Poissonian, the mean of 
the series equals the variance £(/)k(0) = (J . 2 an d T|(/) is 
equal to zero. 

The second step is to specify the “average size” of the 


clusters. The “average” cluster size, defined as the coher- 
ence length, is the length at which the autocovariance 
equals 1/e, where e is the base of the natural logarithm 
and 1/e = 0.3679. The autocovariance is defined as 



where o 2 is the variance for the data series. 

The first step determines the presence of clustering, 
and the second step determines the coherence length or 
the “average size” of the clusters. In the final step, a clus- 
tering intensity parameter is determined as suggested by 
Jameson and Kostinski 4 , defined as: 




(3) 


The clustering intensity parameter depends on the 
value of the two-point autocorrelation function at zero 
lag length and the mean and variance associated with the 
data series. The clustering intensity provides insight into 
the magnitude of the values associated with clustering 
relative to the mean (see ref. 4, Figure 4). In a Poisson 
distribution, where clustering does not occur, the cluster- 
ing intensity is zero. This is because, in a Poisson distri- 
bution, the mean equals the variance, so the last term in 
equation 3 will equal one, making the corresponding clus- 
tering intensity zero. The clustering intensity increases as 
the variance increases relative to the mean for a given 
value of the two-point autocorrelation function. Two data 
series can have the same mean, but the series that con- 
tains values that differ the most from the mean will have 
the greatest clustering intensity. 
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Inflight Measurements 

Liquid water measurements were made by the NASA 
Glenn Research Center’s Twin Otter research aircraft 
during the Supercooled Large Droplet Research Program 
(SLDRP) 14 , conducted from the winter of 1996-97 to De- 
cember 1999. The flight program was a joint effort between 
NASA, the National Center for Atmospheric Research (NCAR), 
and the FAA. The primary objective of the program was to 
acquire flight research data on supercooled large droplets 
(SLDs, defined as larger than 50 pm in diameter) for char- 
acterization of the SLD environment aloft. 

Based on weather forecasts and real-time inflight guidance 
provided by NCAR, the NASA Glenn Icing Research Aircraft 
was flown from the NASA Glenn Research Center to locations 
where conditions were believed to be conducive to the forma- 
tion of SLDs aloft. Research flights were conducted over 
regions adjacent to Lake Erie and Lake Michigan, ex- 
tending as far south as Parkersburg, West Virginia. 
Onboard instrumentation was used to record meteorologi- 
cal characteristics encountered during the flight 14 . Though 
SLD conditions were often located, icing conditions with 
cloud-size droplets less than 50 pm in diameter and mixed 
phase conditions were encountered as well. This occurred 
often when in transit to or from an area with forecasted 
SLDs. 

Liquid water content values used in this study were 
measured by the Twin Otter with a CSIRO-King hot-wire 
probe mounted on the fuselage upstream of the windshield, 
with a range of 0-1 .0 g m \ The probe is sampled at 10 
Hz, and liquid water content is calculated for 1-s incre- 
ments. Probe drift (zero offset) was adjusted by compar- 
ing the King probe mean LWC during a flight segment 
against a Forward Scattering Spectrometer Probe (FSSP) 
mean LWC for the same period. The King probe zero- 
offset was adjusted to allow the King and FSSP mean 
LWCs to match as closely as possible while maintaining 


non-negative LWC values from the King probe. Probe 
zero-offset affects the coherence length and the cluster- 
ing intensity parameter because of the relationship with 
the mean LWC (eq 1-3). 

NASA provided CRREL with LWC measurements for 
32 flights during the winters of 1996-97 and 1997-98. 
Dependence of LWC variability on height above the cloud 
base (or below the cloud top), temperature, presence or 
absence of ice crystals, and turbulence 15 led us to seek 
flight segments that, as much as possible, occurred at a 
nearly constant altitude. In addition, we sought segments 
that were along a nearly constant heading, though we 
often accepted flights of changing direction if they 
did not double back on themselves. For this study 
we chose to ignore the effects of temperature and 
ice crystal variation, which could partially control clus- 
tering. We had insufficient information to exclude turbu- 
lence effects. 

All 32 NASA flights were plotted to show the loca- 
tion of LWC with regard to flight altitude and heading. 
We chose 12 flight segments for analysis because they 
had no breaks in liquid water content and were at nearly 
constant elevation and heading. We also chose long seg- 
ments to represent clustering conditions as faithfully as 
possible. The characteristics of each selected flight seg- 
ment are listed in Table 1 . Flight segment durations ranged 
from 17.4 minutes (73 km) to 45.3 minutes (190 km), 
and mean LWCs ranged from 0.03 to 0.35 g nr 3 . 

Analysis 

As discussed in the methodology, we used a graphical 
threshold approach to show the frequency and percent- 
age of time in each flight segment when the LWC ex- 
ceeds specific threshold values. Figure 2, a hypothetical 
random series with no clustering, showed the time with 
values greater than the thresholds and the number of cross- 




Table 1. Flight segment characteristics. 


Flight 

Duration 

LWC mean 

LWC variance 

Coherence length 

Cluster 


(s) 

(g rrr 3 ) 

(g nr 3 ) 

(s) 

parameter 

97-01 -24f1b 

1495 

0.05 

0.0006 

55 

0.24 

97-01 -24f2a 

1141 

0.08 

0.014 

80 

2.42 

97-01 -24f 1c 

1044 

0.14 

0.008 

147 

0.41 

97-03-14 

2520 

0.16 

0.34 

247 

1.29 

98-02-04a 

2536 

0.03 

0.00006 

8 

0.06 

98-02-04b 

1353 

0.18 

0.0037 

8 

0.11 

98-02-04C 

1156 

0.21 

0.0123 

20 

0.31 

98-02-04d 

1151 

0.19 

0.008 

30 

0.24 

98-02-05 

1709 

0.09 

0.004 

57 

0.42 

98-02-12 

2719 

0.35 

0.10 

225 

0.08 

98-02-24a 

1161 

0.20 

0.015 

149 

0.36 

98-02-24b 

1809 

0.23 

0.008 

25 

0.15 
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ings as being linear and symmetrical, respectively, with 
threshold magnitude. This is not the case with actual cloud 
measurements. The cumulative distance for all cases is 
nearly log-linear, rather than linear, and the number of 
crossings varies considerably with threshold value 
(Fig. 3). 

Comparison of the 12 flight segment plots (not all 
shown here) showed no consistent relationship between 
cluster parameter magnitude and either frequency or per- 
centage of time when LWC is greater than a threshold. 
The plots do suggest, however, that there may be a rela- 
tionship between coherence length and the shape of the 
curve representing the number of threshold crossings with 
threshold magnitude. Flight segments with coherence 
lengths of less than 147 s (10.3 km) all had unimodal dis- 


tributions (Fig. 3a), whereas flight segments with coher- 
ence lengths of 147 s or greater showed multi-modal dis- 
tributions (Fig. 3b). In addition, the plots show evidence 
of physical consistency within given meteorological syn- 
optic situations. For example, flight segments 98-02-04a- 
d all occurred during one flight, and the relationships be- 
tween the threshold and the frequency and percentage of 
time that conditions exceed thresholds are very similar 
for the four segments (see Figure 4 for three of the four 
flight segments). However, this was not true for the mul- 
tiple segments of flights 98-02-24 and 97-01-24. Though 
the cluster parameter was low for all segments of flights 
98-02-04 and 98-02-24, it varied widely in segments of 
flight 97-01-24. We have not analyzed sufficient flights 
to develop a catalogue of distinctive conditions. How- 



Figure 3. Flight segment 98-02-05 (top, Figure 3a) with a coher- 
ence length of 57 s, and flight segment 98-02-24a (bottom, Figure 
3b) with a coherence length of 149 s. 
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Liquid Water Thresholds (g m ' 3 ) 



Liquid Water Thresholds (g m ' 3 ) 

Figure 4. Flight segments 98-02-04b-d (top to bottom, respec- 
tively), all from the same flight, showing remarkably consistent 
crossing and cumulative distance patterns with threshold. 
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ever, we believe that certain characteristics of the vari- 
ability of the LWC along the flight path may be system- 
atically determined with plots of this type. 

The analytical technique indicates that the Jameson 
and Kostinski 4 cluster parameter varies widely, from a 
low of 0.06 indicating almost no clustering, to a high of 
2.42 (Table 1). Jameson and Kostinski 4 provide an indi- 
cation of cluster appearance with cluster parameter mag- 


nitude in their Figure 4. Our Figure 5 shows 1000 s of 
four flight segments covering the range of cluster magni- 
tudes computed for our 12 flight segments. They com- 
pare well to Jameson and Kostinski’s example, and they 
indicate that there is a relationship between cluster pa- 
rameter and the graphical appearance of clustering. 

The coherence length, as shown in Table 1, is mea- 
sured as time in seconds and represents the typical size of 







Figure 5. From top to bottom, flight segments 98-02-04a, 98-02-05, 97- 
03-14, and 97-01-24f2a, showing increasing visual clustering as clus- 
ter parameter increases. The x and y scales on all plots are the same 
for comparative purposes. 
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Figure 6. Demonstration of coherence length computation along 
autocovariance curves for two flight segments. 


clusters. Figure 6 illustrates how the coherence length was 
computed from the He distance along the autocovariance 
curve for two flights. Considering the average Twin Otter 
flight speed of about 70 m s , clusters during the 12 flight 
segments ranged from about 560 m to 17.3 km in length. 
Figure 5 also suggests that there is a relationship between 
coherence length and the visual appearance of cluster size. 

Discussion 

Jameson and Kostinski 4 analyzed LWC in icing con- 
ditions using a simulated series of cloud measurements. 
Based on our small sample size of 12 flight segments from 
only 7 of the 32 NASA SLDRP flights, our analyses indi- 
cate that clustering of actual cloud LWC varies over a 
wide range of cluster intensities and coherence lengths. 
Our largest cluster intensity — 2.42 — exceeds the largest 
value illustrated in Jameson and Kostinski’s Figure 4 — 
2.0 — though they do not provide a maximum cluster value 
for natural conditions. 

We found the cluster intensity to vary over a wide range 
among flights and occasionally among discontinuous seg- 
ments within the same flight. The graphical analyses and 
analytical results show considerable consistency between 
segments of the 98-02-04 flight. This suggests that the 
atmospheric dynamics controlling the clustering persisted 
for a considerable time and distance. Only analyses of 
more flight segments will allow us to judge this better. 
Since one of our goals is to determine if cluster 
climatologies, or characterizations, can be developed for 
regions, for time periods, or for meteorological condi- 
tions such as warm front overrunning or cold air advec- 


tion over water, a large sample of measurements will be 
required to adequately characterize the range of cluster 
conditions and their variance. 

NASA SLDRP flights were intended to locate SLDs, 
so all flights were made into conditions where weather 
forecasters assessed SLDs as a possibility. However, since 
large portions of some flights were made in transit, many 
analyses include non-SLD conditions. We have not as- 
sessed whether SLDs were present in the 12 flight seg- 
ments analyzed, nor have we assessed the meteorological 
conditions documented by NASA and NCAR for the 
flights encountering SLDs. The SLD and meteorological 
information may provide further guidance for explaining 
variation in clustering among flight segments. 

Other factors that may affect clustering intensities and 
coherence lengths include the position of the aircraft with 
respect to the height above the cloud base, the direction 
of aircraft movement with regard to any lineation within 
the cloud structure, turbulence, temperature, and fluctua- 
tions in glaciation from cloud to cloud. We will assess 
some of these factors from available information as we 
analyze additional flight segments to attempt to explain 
segment-to- segment variance. 

Though the graphical presentations of frequency and 
duration above thresholds have value in the analyses, we 
have yet to learn how to use them fully. However, the 
clustering intensities developed from Jameson and 
Kostinski’s 4 methodology, as well as coherence lengths, 
have great value in analytically assessing clustering. Me- 
teorological synoptic conditions, region, and season are 
likely to affect cloud clustering characteristics, but the 
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extent is not known. However, we believe that climato- 
logical characterizations of cloud microphysical condi- 
tions can be constructed, using the methodologies dem- 
onstrated in this paper, to represent the clustering of liquid 
water content, drop and particle size distributions, and 
glaciation. These climatologies could then be used to as- 
sess aircraft icing conditions and the potential performance 
of icing remote sensing systems synoptically, regionally, 
and seasonally. 

Inversion of cluster parameter and correlation length, 
for an assumed LWC distribution for the clustered condi- 
tions, would generate a correlated data series for use by 
modelers assessing icing remote sensing system perfor- 
mance and specifications. Sensitivity analyses could be 
performed by varying coherence lengths, and clustering 
intensities, to generate a spectrum of correlated data se- 
ries that reflect natural conditions. These series could be 
used to investigate the impact of clustered conditions on 
aircraft icing and on radar and microwave radiometer sys- 
tems being proposed for remote detection of aircraft ic- 
ing conditions. 

Though our analyses are of conditions found in one of 
the most icing-intensive areas of North America 16 - 17 , the 
NASA measurements represent only that one region. Pre- 
vious programs assessing icing cloud microphysics, such 
as the Winter Icing and Storms Project (WISP), conducted 
by NCAR in the Front Range area of Colorado 18 , may 
show considerably different clustering conditions. It is 
likely that clustering conditions vary significantly from 
region to region, just as other icing cloud parameters do 7 . 
This paper is a demonstration of the problem as it ad- 
dresses remote sensing needs. 

Conclusions 

Cloud cluster analysis is an element of cloud micro- 
physics that has not received much attention. Clouds are 
neither temporally static or uniform nor spatially uniform. 
However, the difficulty is how to measure and quantify 
variation. Methodologies presented here may be useful for 
assessing clustering for our remote sensing applications. 

In addition to the examples presented here, we will 
assess more cases from the NASA SLDRP, from the 
NCAR Winter and Icing Storms Project (WISP), and from 
the Mt. Washington Icing Sensors Project (MWISP) 19 , 
and as resources permit, we will examine other param- 
eters such as drop size, temperature, particle concentra- 
tion, and glaciation. It may ultimately be useful to assess 


vertical fluctuations of microphysical conditions within 
clouds, a problem more difficult than horizontal varia- 
tion because of the difficulty of acquiring vertical mea- 
surements of cloud microphysical properties and the po- 
tential for greater variation than observed in the horizontal. 

It is likely that most icing clouds, especially those with 
more dangerous conditions, are typically “patchy” and 
rarely Poissonian. Clustering of cloud microphysical pa- 
rameters affects icing processes and potentially the abil- 
ity of remote sensing systems to identify icing conditions. 
Clustering quantification, and perhaps creation of repre- 
sentative data series inverted from climatological sum- 
maries of clustering, will allow more accurate modeling 
of icing and remote sensing systems. CRREL, NASA, 
and the FAA are quantifying clustering of NAS A SLDRP 
flights and developing inversion techniques for these pur- 
poses. 
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